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Understanding  the  dynamics  of  the  glass  transition  of 
polymers  and  molecular  liquids  has  been  a  goal  of 
physicists  for  more  than  a  half-century.  The  nature  of  this 
transition  is  directly  linked  to  the  heterogeneity  of  the 
dynamics,  reflected  in  spatial  and  temporal  correlations  of 
molecular  motions,  as  well  as  a  distribution  of  relaxation 
times.1  4  The  former  can  be  quantified  by  a  dynamic 
correlation  length  that  grows  on  approaching  the  glassy  state 
and  the  latter  from  the  breadth  of  the  relaxation  function  or 
dynamic  susceptibility.  The  majority  of  experimental  studies  of 
this  topic  rely  on  dielectric  spectroscopy,  which  provides 
relaxation  spectra  over  many  decades  of  frequency.  However, 
the  higher  order  susceptibilities  involving  both  spatial  and 
temporal  correlations  cannot  be  obtained  from  linear  relaxation 
measurements  without  assumptions.5  Characterizing  the  spatial 
distribution  of  the  dynamics  requires  other  methods. 

Since  the  1990s  quasi-elastic  neutron  scattering  (QENS)  has 
been  widely  employed  to  obtain  space-specific  information  on 
the  dynamic  properties,  providing  momentum  ( Q)  and 
frequency  dependences  from  the  double  differential  scattering 
function  of  protons.6  The  Q-dependent  scattering  from  glass¬ 
forming  materials  can  be  described  using  the  Kohlrausch— 
Williams— Watts  (KWW)  equation7  11 

^kwwCQ.)  0  K  exp (~[t/Ta(Q_)]m))  (l) 

For  polymers  Ta  is  the  local  segmental  relaxation  time, 
describing  the  time  scale  of  correlated  conformational 
transitions  of  backbone  bonds.  Note  that  the  narrow  dynamic 
range  of  neutron  scattering  introduces  uncertainty  in  the  value 
determined  for  the  stretch  exponent  /?.  For  small  Q,  in  the 
range  below  the  Q  value  (~1  A-1)  for  which  the  static  structure 
factor  shows  its  first  maximum,  tJ1  ~  Q-2, 12,13  the  microscopic 
motions  follow  simple  diffusion  (“Gaussian  approximation”).  At 
larger  Q_  deviations  from  Gaussian  behavior  transpire,  with  the 
weaker  wavevector  dependence  interpreted  by  a  variety  of 
models,  including  mode-coupling  theory  (MCT),  jump 
diffusion  models,  and  caging— decaging  fractals.12,14,15  Molec¬ 
ular  dynamics  simulations  have  also  proven  useful  in 
interpreting  QENS  experimental  data.16 

In  this  work  we  study  the  segmental  dynamics  of  1,4- 
polybutadiene  (PBD)  in  the  presence  of  a  high  concentration 
of  reinforcing  filler  particles.  Reinforcement  of  rubbery 
polymers  is  especially  intriguing  because  hard  particles  increase 
both  the  modulus  and  the  failure  strain  (unlike  plastics,  which 
are  made  more  brittle).  Many  applications  of  elastomers  would 
not  exist  without  the  use  of  particulate  reinforcement  of  the 
material.  The  PBD  was  synthesized  by  anionic  polymerization 
(weight-average  molecular  weight  =  158  kDa;  polydispersity  = 
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1.07;  14%  vinyl  repeat  units;  calorimetric  Tg  =  180  K).  The 
carbon  black  was  N110  (Cabot  Vulcan  9),  comprised  of  17  nm 
particles  fused  into  aggregates  having  a  mean  diameter  of  ca.  54 
nm.  The  aggregates  have  a  highly  structured  surface,  with  an 
effective  volume  much  larger  than  the  actual  volume.  Carbon 
black  has  a  heterogeneous  topography,  with  local  surface 
energies  that  can  be  as  high  as  100  mj/m2.17  N100  is  a  furnace 
black  with  reactive  double  bonds  at  the  particle  surface  that 
serve  as  binding  sites  for  polymer  chains.  After  mixing  on  a 
two-roll  mill,  the  PBD  composite  samples  were  Soxhlet 
extracted  with  cyclohexane  to  remove  unattached  polymer. 
The  initial  composition  of  20.5%  by  volume  carbon  increased 
to  60%  on  extraction  (compositions  determined  by  thermo¬ 
gravimetry).  From  the  surface  area  of  the  carbon  particles  and 
their  size  and  density,18  we  calculate  a  specific  surface  area  of 
4.57  X  105  cm2  per  g  of  composite  (assuming  no  agglomeration 
during  the  extraction);  this  is  about  twice  the  specific  volume  of 
the  PBD  divided  by  its  radius  of  gyration  (~14.6  nm19,20). 
Thus,  for  a  single  layer  having  a  thickness  equal  to  twice  the 
radius  of  gyration,  we  estimate  that  only  about  a  quarter  of  the 
particle  surface  is  covered;  moreover,  every  chain  is  adsorbed 
onto  a  particle  since  the  material  was  extracted. 

Interaction  and  binding  of  polymer  segments  to  the  particles 
constrain  their  motions  increasing  fa,21  and  molecular 
dynamics  simulations  have  shown  that  interactions  with 
nanoparticles  increase  the  length  scale  of  the  cooperative 
motions  of  polymers  underlying  ra.22  Interestingly,  QENS  of 
polyisoprene  below  its  glass  transition  temperature  showed 
increased  librational  mobility  of  the  backbone  in  the  presence 
of  carbon  black,”3  although  whether  this  result  is  specific  to  the 
glassy  state  remains  to  be  investigated.  More  generally,  the 
spatial  variation  of  the  constraints  from  particular  fillers  is 
unclear,  and  in  order  to  explore  this  aspect  of  the  dynamics,  we 
employed  quasielastic  neutron  scattering.  Experiments  were 
carried  out  using  the  high  flux  back-scattering  spectrometer 
(HFBS,  NG2)  at  the  NIST  Center  for  Neutron  Research.  The 
dynamic  scattering  function,  S(Q,E),  was  measured  with  an 
energy  resolution  of  0.4  /<eV  at  up  to  17  fieV;  the  Q  range  was 
0.31-1.71  A-1.  A  sample  thickness  of  0.1  mm  was  chosen  to 
avoid  multiple  scattering.  For  our  samples,  we  estimate  the 
incoherent  scattering  cross  section  of  hydrogen  atoms  in  the 
polymers  to  be  94%  of  the  total  scattering  for  neat  PBD  and 
82%  for  the  extracted  nanocomposite.  Therefore,  the  scattering 
spectra  reflect  the  energy-dependent  distribution  of  self¬ 
correlation  of  the  hydrogen  atoms. 
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The  desired  KWW  function  was  expressed  as  a  weighted 
superposition  of  Lorentzian  functions 

$kwwM  =  //(lo g  r)  WT 2  2  d(logr) 

J  1  +  CO  T  (2) 

For  the  distribution  function,  /(log  r),  we  used  the  generalized 
version24  of  the  expression  of  Rajagopal  and  Ngai.25  This  model 
function  was  first  convoluted  with  the  experimental  resolution 
function  and  then  fit  directly  to  the  experimental  S(Q,ft>)  to 
obtain  r„(Q)  and  /?(Q)  at  the  various  measurement  temper¬ 
atures. 

Figure  1  shows  master  curves  of  the  dynamic  mechanical  loss 
modulus,  measured  in  shear  using  an  Anton  Paar  MCR  502 


Frequency  [Hz] 

Figure  1.  Master  curves  of  the  mechanical  loss  moduli  of  the  PBD  neat 
and  with  20.5%  by  volume  carbon  particles  along  with  the  fitted  KWW 
function.  The  inset  shows  the  frequency— temperature  shift  factors, 
which  have  the  same  temperature  dependence.  Note  for  these 
measurements  only,  the  filled  sample  was  not  extracted.  Error  bars  are 
not  larger  than  the  symbol  size. 


rheometer,  for  PBD  and  the  unextracted  nanocomposite.  The 
spectrum  for  the  neat  polymer  is  well-described  by  eq  1  with  /? 
=  0.41,  which  is  in  good  agreement  with  the  value  determined 
by  dielectric  spectroscopy.26  Determining  /?  for  the  filled 
polymer  is  problematic.  The  mechanical  loss  peak  in  Figure  1  is 
significantly  broader  than  for  the  neat  PBD  and  can  be 
approximately  described  using  /?  =  0.38,  with  deviation  at  lower 
frequencies.  The  presence  of  filler  heterogeneously  broadens 
the  segmental  relaxation  because  the  effect  of  the  carbon 
particles  on  the  dynamics  depends  on  the  proximity  of  the 
segments  to  the  particle.  Of  course,  neither  mechanical  nor 
dielectric  spectroscopy  can  provide  information  on  this  In¬ 
dependence;  the  loss  modulus  measurements  average  over  all 
length  scales. 

Figure  2  compares  representative  dynamic  structure  factors 
for  the  two  materials  at  T  =  300  K.  S(E )  for  the  nanocomposite 
is  narrower,  which  indicates  suppression  of  the  dynamics.  Fits 
of  eq  2  to  the  spectra  in  Figure  2  are  shown.  For  the  neat  PBD, 
/?  =  0.41,  with  the  (^-dependence  of  relaxation  times  shown  in 
Figure  3a  for  four  temperatures  above  Tg.  For  Q.  <  1  A-1  the 
Gaussian  approximation  holds,  t„°'41  ~  Q.-2.  Note  that  the 
stretch  exponent  for  the  incoherent  dynamic  structure  factor  is 
(^.-independent  and  thus  is  the  same  as  measured  by 
mechanical  and  dielectric  spectroscopies.  Included  in  Figure  3a 
are  results  from  ref  13  for  Q.  <  1  A-1;  the  agreement  with  the 
present  data  is  satisfactory.  At  shorter  length  scales  Fickian 
diffusion  does  not  apply,  and  the  Q-dependence  becomes 


Figure  2.  Incoherent  dynamic  structure  factors  at  300  K  and  the 
indicated  wavevectors  for  neat  PBD  (circles)  and  the  nanocomposite 
(triangles)  along  with  fits  to  eq  2  (solid  lines).  The  dashed  lines 
indicate  the  experimental  resolution  measured  at  4  K 
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Figure  3.  Local  segmental  relaxation  times  from  QENS  for  (a)  unfilled 
1,4-polybutadiene,  with  the  value  of  the  stretch  exponent  obtained 
from  mechanical  spectroscopy  (see  Figure  l);  (b)  extracted  nano¬ 
composite  with  the  /?(Q)  (plotted  in  the  inset)  taken  to  he  the  values 
giving  conformance  to  the  Gaussian  approximation  for  <0.8  A-1  and 
fixed  at  0.41  for  higher  Q.  Lines  are  least-squares  fits  having  a  slope 
equal  to  —2  at  small  Q. 


weaker,  Ta0'41  ~  QT°'9,  indicating  more  local  and/or  isolated 
motions.  The  wavevector  associated  with  this  crossover  in 
behavior  is  the  same  for  all  temperatures. 

As  noted,  we  cannot  obtain  from  mechanical  spectroscopy 
the  value  of  /?(Q)  for  the  nanocomposite;  instead,  we 
determine  at  each  Q_  the  f)  value  yielding  the  expected  za  ~ 
q-2//'  Behavior  for  small  Q.  (recognizing  that  the  Q.  dependence 
of  ta  may  depend  weakly  on  /?).  These  results  are  shown  in 
Figure  3b;  note  that  the  f)  are  independent  of  temperature. 
Similar  to  the  results  for  the  neat  PBD,  there  is  a  transition  to  a 
very  weak  Q_  dependence  for  larger  Q,  although  this  crossover 
occurs  at  a  somewhat  smaller  wavevector  (larger  spatial  range) 
than  for  neat  PBD.  The  /?(Q)  for  the  nanocomposite  are 
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plotted  versus  Q.  in  the  inset  to  Figure  3b.  For  larger  Q,  all 
segments  are  constrained  by  the  particles,  and  the  distribution 
of  relaxation  times  is  narrower  than  at  larger  length  scales.  The 
increase  in  the  mean  relaxation  time  due  to  the  filler  particles  is 
larger  at  larger  Q.  For  example,  for  the  data  in  Figure  2,  at  Q  = 
0.68  A-1  ra  is  about  6-fold  larger  for  the  nanocomposite 
compared  to  the  neat  PBD,  whereas  at  Q.  =  1.56  A-1  the 
segmental  relaxation  increases  a  factor  of  36  due  to  the  filler 
particles. 

In  conclusion,  binding  of  PBD  chains  to  the  carbon  particle 
surface  slows  down  the  segmental  motions  and  increases  the 
heterogeneity  of  the  dynamics,  with  the  spectral  stretching  and 
change  in  ra  becoming  spatially  dependent;  that  is,  the  effect  of 
the  particles  depends  on  the  spatial  distance  of  the  segments, 
notwithstanding  that  all  chains  are  attached  to  the  particles. 
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